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Abstract 

We investigate observational constraints on the normal branch of the warped DGP 
braneworld cosmology by using observational data from Type la Supernovae (SNIa), 
Baryon Acoustic Oscillations (BAO), Cosmic Microwave Background (CMB) and Baryon 
Gas Mass Fraction of cluster of galaxies. The best fit values of model free parameters are: 
Vtm = 0.240^^30 and = 0.000+°°^^ at la confidence interval by using Gold sample 
SNIa+CMB shift parameter+BAO+Gas mass fraction of baryons in cluster of galaxies. 
The results for essence sample SNIa combined with CMB shift parameter, BAO and 
Baryon Gas mass fraction correspont to: = 0.220+|];?7[J and n^.^ = 0.000+°-°25 at 68.3% 
confidence interval. We determine the age of the universe by using these best fit values. 
We also study the effective cosmological dynamics on the brane via an effective equation of 
state parameter and the deceleration parameter to conclude that an effective phantom-like 
behavior arises in this scenario. 
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1. Introduction 

The accelerated expansion of the universe 
supported by recent observational data [1] 
is one of the most important discoveries in 
the last decade for the cosmology commu- 
nity. Within the framework of the gen- 
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eral relativity, the acceleration could be 
associated with the so-called dark energy, 
whose theoretical nature and origin are still 
unknown for theorists. Cosmological con- 
stant or vacuum energy with an equation 
of state parameter u = —1, is the most 
popular candidate for dark energy but un- 
fortunately, it suffers from some serious 
problems such as fine-tuning and coinci- 
dence problems. Therefore, a number of 
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models containing dynamical dark energy 
have been proposed as the mechanism for 
late-time cosmic speed up [2]. Some of 
them are quintessence, k-essence, phantom 
scalar field, chaplygin gas models and so 
on. Another alternative approach to explain 
the late-time cosmic speed up is modifica- 
tion of the geometric sector of the Einstein 
field equations leading to modified gravity 
[3]. In the spirit of modified gravity pro- 
posal, the Dvali-Gabadadze-Porrati (DGP) 
braneworld scenario explains the late-time 
accelerated expansion in its self-accelerating 
branch without need to introduce a dark 
energy component on the brane. In this 
scenario our universe is a 3-brane embed- 
ded in a higher dimensional flat space-time 
(bulk). The late-time acceleration is driven 
on large scales by leakage of gravity from 
the brane into the bulk ( see for instance 
[4] and references therein) . On small scales, 
gravity is bound to the brane and general 
relativity is recovered to a good approxima- 
tion. Indeed, the DGP model has two types 
of solutions (corresponding to two possi- 
ble embedding of brane in the bulk ): the 
self-accelerating, (DGP^'^^) branch and the 
normal, {DGP^-^) branch. The (DGP^^^) 
branch however, suffers from some insta- 
bilities such as ghosts [5] and cannot de- 
scribe the early stages of the universe evolu- 
tion properly. There are some extension of 
the DGP setup that provide relatively wider 
parameter spaces with richer phenomenol- 
ogy. One of these models is the warped 
DGP braneworld (WDGP) [6] which is a 
unified model of Randall- Sundrum II (RSII) 
braneworld scenario [7] and the DGP setup. 
In the RSII braneworld model, gravity mod- 
ifies in early ( high energy) epoches of the 
universe evolution. The warped DGP sce- 
nario also gives a self-accelerating phase in 
the brane cosmology. It is important to note 



that, the self-accelerating branch gives an 
effective equation of state that never can be 
less than —1 (always non-phantom behav- 
ior and therefore no crossing of the phan- 
tom divide line). The other branch of this 
scenario is the normal branch {DGP^~^) 
which doesn't self-accelerate but requires 
dark energy or modification of the induced 
gravity on the brane to explain the late 
time acceleration. It is possible to real- 
ize phantom-like effects (without phantom 
matter) in this normal branch via screening 
of the brane cosmological constant at late 
time [8]. 

This paper is devoted to explore obser- 
vational status of the normal branch of the 
warped DGP braneworld cosmology. Ob- 
servational constraints in DGP model with 
and without tension is investigated in Ref. 
[9]. We impose constraints on the model pa- 
rameters by using the several recent obser- 
vations such as distance measurements from 
type la supernovae (SNIa) from the Gold 
[10] and Essence [11] surveys, the baryon 
acoustic oscillations (BAO) measurement 
from the large-scale correlation function of 
the Sloan Digital Sky Survey (SDSS) [12], 
the position of the first peak of the cos- 
mic microwave background (CMB) from 
WMAP7 [13] and the baryon gas mass frac- 
tion of cluster of galaxies [14]. The struc- 
ture of the paper is as follows: in section 
2, we introduce the model and its cosmo- 
logical implications. In sec. 3 we explore 
the effect of the WDGP model on the ge- 
ometrical parameters of the universe. Sec- 
tion 4 includes the observational constraints 
on the model parameters space while in sec- 
tion 5 the detailed results are presented. In 
section 5 we perform a detailed comparison 
between age of the oldest objects of the uni- 
verse and the age result that obtained from 
the best fit values of our model parameter 



space. We also study the effective dynamics 
of the model in this section. Finally section 
6 is devoted to the concluding remarks. 



2. The Model 

The action of 
brane world model 
follows [6,15] 



the warped DGP 
can be written 
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Here Sbuik is the action of the bulk, Sbrane 
is the action of the brane and S is the total 
action. with A = 0,1, 2, 3, 5 are coordi- 
nates in the bulk, while x'^ with /i = 0, 1, 2, 3 
are induced coordinates on the brane. k| is 
5-dimensional gravitational constant. ^^^R 
and ^^^Cjri are 5-dimensional Ricci scalar 
and matter Lagrangian respectively. 
is trace of the extrinsic curvature on either 
sides of the brane. Cbraneidapyip) is the ef- 
fective 4-dimensional Lagrangian. The ac- 
tion S is actually a combination of the 
Randall- Sundrum II and the DGP model. 
In other words, an induced curvature term 
is appeared on the brane in the Randall- 
Sundrum II model. Now we consider the 
brane Lagrangian as follows 

C-braneidap-, = — -R — A + L^, (3) 

where /i is a mass parameter, R is the Ricci 
scalar of the brane, A is the tension of the 
brane and is the Lagrangian of the other 
matter fields localized on the brane. We as- 
sume that bulk contains only a negative cos- 
mological constant, A5. With these choices. 



action (1) gives either a generalized DGP 
or a generalized RS II model: it gives DGP 
model if A = and A5 = 0, and gives RS II 
model if = 0. The generalized Friedmann 
equation on the brane is as follows [6] 
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where e = ±1 is corresponding to two possi- 
ble branches of solutions (two different em- 
bedding of the brane) in this warped DGP 
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m\ + 
ms = 



k^'^^^ . Also, Sq is an integration constant 
and corresponding term in the generalized 
Friedmann equation is called dark radiation 
term. We neglect dark radiation term in 
what follows. In this case, the generalized 
Friedmann equation (4) takes the following 
form 
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where p = pm is the energy density of dark 
matter on the brane. 

2.1. Cosmological Implications 

In this section we study cosmological dy- 
namics on the DGP brane embedded in a 
warped bulk Manifold. To this end, we as- 
sume a flat FRW universe on the warped 
DGP brane. In this setup we can rewrite 
equation (5) as follows 




where 



\ is the DGP crossover scale. 



2m^ 

In the distance scale lower than this scale, 
gravity behaves as usual general relativis- 
tic one but in the distance scales higher 
than the crossover scale, gravity leaks to 
the extra dimension and this leakage leads 
to weakness of gravity in the large scales, so 
the universe expansion accelerates. The up- 
per sign in equation (6) corresponds to the 
self-accelerating branch of the model. Tak- 
ing the lower sign of this equation results a 
very interesting feature. Indeed, if we as- 
sume a model universe with standard cold 
dark matter (SCDM) with pm = Po„ (^) 
the accelerating behavior of the model can 
be recovered by rewriting Friedmann equa- 
tion (6) as follows 



+ A, 



where Ag// mimics the role of an effective 
cosmological constant on the brane (note 
that it is not actually a constant!) and it 
can be decomposed into two parts as follows 
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The first two terms appeared in parenthe- 
sis on the right hand side of this relation 
could be considered collectively as a cosmo- 
logical constant term on the brane and the 
last term on the right hand side (the square 
root) screens the effect of the brane cos- 
mological constant in the same way as has 
been pointed out by Lue and Starkman [8]. 
In this situation, the effective cosmological 
constant A^ff on the brane increases with 
time due to dynamical screening effect, that 
is, reduction of the second term on the right 



hand side of (8) with cosmic time. In fact 
the normal branch of the model has the key 
property that brane is extrinsically curved 
so that shortcuts through the bulk allow 
gravity to screen the effects of the brane 
energy-momentum contents at Hubble pa- 



rameters H 



The screening effect is 



a result of leakage of gravity to the extra 
dimension at late times. 
For future purposes, it is useful to express 
the Friedmann equation (6) in a dimension- 
less form as follows 
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lOOh km/s/Mpc. Taking z = imposes 
a constraints on the model parameters as 
follows 
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The general relativistic limit can be recov- 
ered if we set = (or = 0). In this 
case equations (10) implies that Q^ + ^x = 
1. Considering the tension of the brane 
cosmological constant this case leads to the 
ACDM cosmology. 

3. The effect of WDGP on the Ge- 
ometrical Parameters of the Uni- 
verse 

The cosmological observations are mainly 
dependent on the background geometry ( 
especially background spatial curvature) of 
the universe. So, in this section we study 
the effect of the WDGP model on the 
geometrical parameters of the universe. 



1. Comoving radial Distance 

One of the basic parameters in cosmology 
is the comoving radial distance. For an ob- 
ject with redshift 2; in a FRW background, 
this parameter can be expressed as follows 
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where J-" = {x, sinx, sinhx) for K = 
(0, 1, —1) respectively, i^' marks curvature 
of the spatial geometry and = '^jjjp- 
Figure 1 shows the radial comoving distance 
versus the redshift for different values of Vl^^ 
in a flat background. Clearly, increasing the 
values of Vtr^ results in a longer comoving 
distance. The mentioned quantity is a use- 
ful quantity in the analysis of the luminosity 
distances of Supernova type la. 
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Figure 1: Comoving radial distance versus the 
redshift for different values of & FRW 

universe. The unit of vertical axis is c/Hq. 

2. Angular size 

The apparent angular size of an object 
located at the cosmological distance is an- 
other important parameter that can be af- 



fected by the cosmological model during the 
history of the universe. If we take the ob- 
ject to lie perpendicular to the line of sight 
and to have physical extent the apparent 
angular size 9 is given by 
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where dA^z) = r{z)/{l + z) is the angular 
diameter distance which is a measure of how 
large objects appear to be in the universe. 
A key application of equation (12) is in the 
study of features of the cosmic microwave 
background radiation. The variation of ap- 
parent angular size AO in terms of the Az 
is given by 



% = H{z)r{z) 



(13) 



This relation is the so-called Alcock- 
Paczynski test. The advantage of the 
Alcock-Paczynski test is that in this case, 
instead of using a standard candle, we 
use a standard ruler such as the baryonic 
acoustic oscillation. Figure 2 shows Az/A6 
for different values of in a flat FRW 
background. As the figure shows, increas- 
ing of results in larger values of Az/ AO. 

3. Comoving Volume Element 

The comoving volume element is another 
geometrical parameter which is used in 
number-count tests such as lensed quasars, 
galaxies, or clusters of galaxies. The co- 
moving volume element can be expressed in 
terms of comoving distance and Hubble pa- 
rameter as follows 
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Figure 3 shows the comoving volume ele- 
ment versus the redshift for different values 



Figure 2: lS.z/ /^O versus the redshift for differ- 
ent values of in a flat FRW background. 

of VLr^ in a flat FRW background. This fig- 
ure indicates that the value of the comoving 
volume element increases with the increas- 
ing of fire- We note the quantities displayed 
in figures 1-3 are specified just by giving 
the value of fi^^ alone. This is because we 
are interested in the DGP character of the 
model. Nevertheless, in plotting these fig- 
ures we have used the values VLm = 0.24, 
Qf^^ = 1 and the other quantity Q\ is ob- 
tained via constraint equation, (10). 

4. Observational Constraints 

In this section, we constrain the model 
parameters of the warped DGP scenario by 
using the most recent observational data 
including the SNIa data measurements as 
given by the gold and essence samples, 
combined with the information from the 
BAO measurements by SDSS, the CMB 
shift parameter from WMAP7 observations 
and the baryon gas mass fraction. 

A: SNIa data 



Figure 3: Comoving volume element versus the 
redshift for different values of in ^ Acit FRW 
background. 

Here we use some new sample supernova 
type la data such as the gold sample com- 
piled in [10] and essence sample compiled in 
[11] to constrain the free parameters of the 
model. This observation directly measures 
the apparent magnitude m of a supernova 
versus its redshift z. We note that deference 
between SNIa samples is in their systematic 
errors and the range of redshift which the 
apparent magnitude is determined . 

The theoretical distance modulus is de- 
fined as 

fi^m-M ^5 logio dL + 5 logio (j~^^f^ + 

(15) 

where M is the absolute magnitude that is 
believed to be constant for all Type la su- 
pernovae and di is the luminosity distance 
that in general case can be expressed as fol- 
lows 

y/Wli] \ Jo E{z'-{ei})J 

(16) 

where {6i\ denotes the model parameters 
and E{z) is given by equation (9). We es- 
timate the best fit of the set of parameters 



{9i} by using a statistics, with 
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In this relation, is the number of SNIa 
data points which is different in several data 
sets, ^obs is the observed distance modulus 
and the aj is the uncertainty in the observed 
distance modulus, which is assumed to be 
Gaussian and uncorrelated so that the like- 
lihood is proportional to exp(— x^/2). The 
parameter M = Slogj^Q d ^^^jpc ) + 25 is a 
nuisance parameter and it is independent of 
the data and the data sets. Following the 
techniques described in Ref. [16], the mini- 
mization with respect to M can be made by 
expanding the Xsn equation (17) with re- 
spect to M as 

xlNiiOi}) = A- 2MB + M^C (18) 
where 

[Mobs(2j; {e^}) - {ej, A/ = o)]2 
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Equation (18) has a minimum for M = 

§ at xiAm = Am) - 

Using this equation the best fit values of 
model parameters as the values that mini- 
mize XsNiaii^i}) can be obtained. 

For the likelihood analysis we marginal- 
ize the likelihood function L ~ exp(— x^/2) 
over h. We adopted Gaussian priors such 
that h = 0.705 from the WMAP7 [13]. 
Table 1 summarizes these priors. 



Table 1: Priors on the parameter space used in the 
Ukehhood analysis. 



Parameter 


Prior 




0.000 - 1.000 Top Hat 




0.000 - 1.000 Top Hat 




-1.000 - 1.000 Top Hat 


h 





(19) 



B: CMB shift parameter 

We use the CMB data from WMAP7 ob- 
servation that includes the shift parameter 
TZ and the redshift of the decoupling epoch 
. The shift parameter TZ relates the angu- 
lar diameter distance to the last scattering 
surface, the comoving size of the sound hori- 
zon at = 1091.3 and the angular scale of 
the first acoustic peak in the CMB power 
spectrum of the temperature fluctuations. 
The CMB shift parameter is approximated 
by (see for more details Ref. [16]) 



7^ 




(21) 



where dA^z) is the angular diameter dis- 
tance defined by equation (12). The con- 
straints on a typical model using CMB shift 
is obtained from minimization of the quan- 
tity 



XCMB 



(22) 

'CMB 

where TZots = 1-725 is the observed value of 
the CMB shift parameter performed from 
WMAP7 observation and its corresponding 
la error is <Jcmb = 0.018 [13]. Also TZth 
corresponds to the theoretical value of shift 
parameter calculated from equation (21). 

C: BAO observation 
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The baryonic acoustic oscillation (BAO) 
peak detected in the SDSS luminous red 
Galaxies (LRG) is another tool to test the 
model against observational data. BAO are 
described in terms of a dimensionless pa- 
rameter 



■Aizsds3;{di}) = vf^ 



Hodl{zsdss;{9i}) 



H{Zsdsa; {^i})2sdss(l + Zsdas) 
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The fo^^ the BAO is given by 
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-0.35 



The observed value from the LRG is 

± 0.017 measured 

at Zsdss = 0.35 [20]. Here Ug = 0.963 is 
the spectral index as measured by WMAP 
seven years observations [13]. 

D: Gas mass fraction of cluster of 
galaxies 

Another cosmological test to constrain 
the parameters of the model arises from 
baryon gas mass fraction of cluster of galax- 
ies for a range of redshifts 
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The basic assumption underlying this 
method is that the baryon gas mass frac- 
tion in clusters is constant, independent of 
redshift. This method can give a constraint 
to the geometry of the universe with the re- 

3 

lation Sgas oc d\ under the assumption that 
this fraction should be approximately con- 
stant with redshift. Following [17] (see also 
[16]), the expression for gas mass fraction 
is given by 



lsfaU^j;{e,})~sl'u^j;{e,})]^ 



(26) 



where Sgas is a dimensionless parameter de- 
fined as 



S, 
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In this relation dJ)^^ is the angular diame- 
ter distance to a cluster in the test model 
which is assumed to be SCDM (cold dark 
matter) in this case, and 6 is a bias fac- 
tor suggesting that the baryon fraction in 
clusters is slightly lower than for the uni- 
verse as a whole. Also 1 + /3 is a factor 
taking into account the fact that the total 
baryonic mass in clusters consists of both 
X-ray gas and optically luminous baryonic 
mass (stars), the latter being proportional 
to the former with proportionality constant 
/3 ^ 0.19^/5^ [16-18]. The nuisance param- 
eter f = TTaTT^ should bc marginalized via 
expanding the x^as equation (26) with re- 
spect to ^ which gives 
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We use the 26 cluster data [14] for 5°^^ to 
obtain the best fit parameters of the model. 

It is important to note that the above 
observational data are uncorrelated since 
they are given by different experiments and 
methods. So, we can construct a joint anal- 
ysis as 



Xt 



ot 



XSN + XCMB + XSDSS + Xlas ■ (31) 



5. Results 



With these prehminaries, we have ob- 
tained the best fit parameters of the normal 
branch of the WDGP model for SNla data 
(gold and essence datasets), the joint anal- 
ysis of the SNIa and CMB, the combined 
analysis of the SNIa, CMB and SDSS and 
finally the joint analysis of the total datasets 
including cluster galaxies gas mass fraction. 
Table 2 shows the results of the observa- 
tional constraints on the free parameters of 
this model. Figure 4 shows the marginal- 
ized relative likelihood with respect to pa- 
rameter and fitted with SNIa gold 
sample, SNIa+CM]3, SNIa+CMB+SDSS 
experiments and SNIa+CMB+SDSS-|-gas 
mass fraction observations. In figure 5, we 
repeat these stages for Essence sample of 
SNIa experiments together with other data 
sets. We plot contour maps of ver- 
sus Qm in figure 6 for Gold Sample and 
combined observational data sets. In fig- 
ure 7 we plot corresponding contour maps 
for Essence Sample and combined observa- 
tional data sets. 
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Table 2: The fitting results for WDGP model by using the SNIa (Gold and Essence Samples), SNIa+CMB 
and SNIa+CMB+SDSS and SNIa+CMB+SDSS+glaxies clusters gas mass fraction experiments in a flat 
background. 



Observation 










Xiiin/^d.o.f 


Agc(Gyr) 


SNIa( Gold Sample) 




u.iou_o.i5o 




0.009 


0.923 


64.97±^^:^^1! 


SNIa(Gold)+CMB 


0.280±U:\'^lj 




0.631±U:f/^, 


0.939 


0.943 


14.38+-- 


SNIa(Gold)+CMB+SDSS 


260+" '^'^" 

U.ZOU_o.l20 


0.000+°°i2 




0.839 


0.992 


i4.5i±^:ti 


SNIa(Gold)+CMB+SDSS+GAS 


U.Z'±U_o.l30 


000+0014 




0.999 


0.992 


14 «f;+i+t)i 

i4.00_3 942 


SNIa( Essence Sample) 


U.UZU_o.015 


170+0.020 


lfi1+0-0V4 


0.009 


1.032 


78 QO+56.018 
(0.00_42 390 


SNIa(Essence)+CMB 


2-50+"-"^" 

U.ZOU_Q J^gQ 


000+0019 


n fiQI +O il'i' 


0.539 


1.030 


14.86+--^ 


SNIa(Essence)+CMB+SDSS 


220+" "^'^ 


00 +0 022 


fiq7+".i4b 
u.Dy(_o.095 


0.689 


1.054 


15.071-- 


SNIa(Essence)+CMB+SDSS+GAS 


220+'^ "^" 

^■^^"-0.170 


0.000 +0 025 


710+0.187 


0.119 


1.044 


14.96±^:i^ 




Figure 4: Marginalized relative likelihood with respect to parameter (left) and fi^c 

(right) fitted with SNIa Gold Sample, SNIa (Gold)+CMB, SNIa(Gold)+CMB+SDSS and 
SNIa(Gold)+CMB+SDSS+ galaxies clusters gas mass fraction experiments. The intersection of 
each curve with horizontal solid line is corresponding to the bound with la confidence level. 
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Figure 5: Marginalized relative likelihood with respect to parameter ilm (left) and ^Ir^ (right) 
fitted with SNIa Essence sample, SNIa (Essence)+CMB, SNIa (Essence)+CMB+SDSS and SNIa 
(Essence)+CMB+SDSS+galaxies clusters gas mass fraction experiments. The intersection of each curve 
with horizontal solid line is corresponding to the bound with Icr confidence level. 




0.1 0.2 0.3 0.4 0.1 0.2 0.3 0.4 0.1 0.2 0.3 0.4 0.1 0.2 0.3 0.4 



Q Q Q. £1 

m m m m 

Figure 6: Contour maps of ^Ir^ versus for WDGP model with la (68.3%) and 2a (95.5%) confi- 
dence levels using SNIa (a), SNIa+CMB (b), SNIa+CMB+SDSS (c) and SNIa+CMB+SDSS+Gas (d). We 
used the Gold Sample for SNIa experiment. 
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Figure 7: Contour maps of ilrc versus 0,m for normal branch of WDGP model with la (68.3%) 
and 2a (95.5%) confidence levels using SNIa (a), SNIa+CMB (b), SNIa+CMB+SDSS (c) and 
SNIa+CMB+SDSS+Gas (d). We used the Essence Sample for SNIa experiment.. 
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5.1. The Age of the universe 

The age of the universe for an expanding 
universe in a flat background is given by 

to= I dt= I (32) 



.0 + 

The WMAP collaboration [13], using the 
AC DM model, quotes to = 13.75 ±0.11 
Gyr. In our model H{z) is given by Eq. 
(9). A model- independent limit of to can be 
obtained from the age of the oldest stars, 
tstar- Since the age of the universe should be 
greater than the age of the oldest stars, this 
limit can be used as another probe to con- 
strain our model ( see for instance [19,20]). 
Studies on the old stars [21] suggest an age 
of 13^^ Gyr for the universe. Richer et. 
al. [22] and Hansen et. al. [23] also pro- 
posed an age of 12.7 ± 0.7 Gyr, using the 
white dwarf cooling sequence method. Re- 
cently, Frebel et al. [24] have reported the 
discovery of HE 1523 — 0901, with an age 
tstar = 13.4±0.8(l(T)±1.8(s?/st) Gyr. With 
these statistical and systematic errors, we 
adopt tstar > 12 Gyr for constraining our 
model [21]. Table 2, shows that the age 
of the universe is always greater than tstar 
for all used data sets. However, the SNIa 
data sets (Gold or Essence Samples) alone 
lead to very great to which is not reliable 
and it seems that the pure SNIa test cannot 
give a good constraint on the model param- 
eters space. This feature is evident also in 
our previous relative likelihood figures and 
has its origin in the geometric nature of the 
SNIa data sets. On the other hand the com- 
bined analysis of SNIa with other data sets 
improves the constraints and gives a more 
reasonable results. 

5.2. Effective Dynamics 

It has been shown in our previous work 
[15] that the Warped DGP scenario has a 



phantom-like behavior without need to in- 
troduce a phantom matter neither in the 
bulk nor on the brane. To investigate the 
effective dynamics of the model we study 
the effective density defined by 



Peff = A4 



2rl 



1 + 4r2(^ 



3/i2 



P + A As 
3 
(33) 

and an effective equation of state parameter 



-1 



Peff 



(34) 



So, the effective equation of state of dark 
energy can be expressed as follows 



= -1- 



(35) 

Using the best fit parameters obtained 
here and summarized in table 2, in figure 
8 we plot ujeff versus the redshift for com- 
bined SNIa+CMB+SDSS+galaxies clusters 
gas mass fraction data set for Gold and 
Essence Samples. 

As the figure shows, in low redshifts and 
especially at present, dark energy exhibits a 
phantom-like acceleration with Ueff < — 1 ( 
it is easy to show also that in this model 
H < and Pe// > 0, see [15] for details). 
In future { z < 0) the effective equation of 
state parameter approaches the cosmologi- 
cal constant line Ugff = — 1. It is important 
to note that although Wg// has a phantom- 
like dynamics in this warped DGP scenario, 
ojtot doesn't show such a behavior. The ex- 
pression for Wtot is given by 



^tot 



(36) 
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Figure 8: The effective equation of state 
parameter versus the redshift fitted with 
SNIa+CMB+SDSS+galaxies clusters gas mass 
fraction with Gold and Essence Samples re- 
spectively. 
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Figure 9: Total equation of state pa- 
rameter versus the redshift fitted with 
SNIa+CMB+SDSS-h galaxies clusters Gas 
mass fraction with Gold and Essence Samples 
respectively. 



Indeed, the total equation of state param- 
eter remains Quintessence-like [wtot > —1) 
as figure 9 shows. 

Another dynamical quantity of interest 
is the deceleration parameter which is ex- 
pressed as g = 



HoE^ 



I where 



E 



1 - 



(37) 



This equation implies that for all values 
of 2;, E < 0. As an important result, the 
deceleration parameter could be never less 
than —1. Consequently, there is no super- 
acceleration or big rip singularity in this 
model. Figure 10 shows the variation of the 
deceleration parameter versus the redshift 
with two combined data sets. It is clear that 
the deceleration parameter reduces by red- 
shift towards recent epoch and in the future 



it approaches q = —1. 

6. Summary and Conclusion 

The normal branch of the warped DGP 
brane scenario has several fascinating prop- 
erties. Its cosmological dynamics does not 
suffer from instabilities such as ghost, it can 
exhibit a phantom mimicry and it leads to 
the famous ACDM scenario in certain lim- 
its. In this regard, we have studied cosmo- 
logical constraints imposed on this model 
from observational data such as type la su- 
pernova data from the Gold and Essence 
surveys, the baryon acoustic oscillations 
(BAO) measurement from the Sloan Digital 
Sky Survey (SDSS), the cosmic microwave 
background (CMB) and the baryon gas 
mass fraction of clusters of galaxies. We find 
that the best fit values of model free param- 
eters are constrained to: VL^ = 0.2401q;53o 
and VLj.^ = 0.000+°'°^^ at la confidence inter- 
val by using Gold sample SNIa -|- CMB shift 
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Figure 10: Variation of the deceleration param- 
eter versus the redshift. The parameters values 
used to produce this figure are taken from table 
2. 



via effective equation of state parameter 
and the deceleration parameter by using the 
best fit parameters of the combined analy- 
sis of SNIa+CMB+SDSS+galaxies clusters 
gas mass fraction with Gold and Essence 
Samples. This analysis confirms that the 
WDGP model has a phantom-like behavior 
without need to introduce phantom fields 
neither in the bulk nor on the brane. It is 
important to note that in this scenario the 
effect of the bulk's warped geometry (en- 
coding in A5) is that by increasing the abso- 
lute value of the bulk cosmological constant 
the phantom-like behavior reduces. In other 
words, incorporation of A5 leads to a reduc- 
tion of the effective phantom nature of the 
model in comparison with pure DGP case ( 
see Ref. [15]). 



parameter + BAO + gas mass fraction. The 
same analysis just including essence sample 
SNIa instead of Gold sample correspond to: 
n^n = 0.2201^70 and = 0.000+°-°25 at 
68% confidence interval. Furthermore, us- 
ing these best fit parameters, we obtained 
the age of the universe 14.85l3;942 and 
14.96l7'g8g at Icr confidence interval by us- 
ing SNIa+CMB shift parameter+BAO+gas 
mass fraction for gold sample and essence 
sample of SNIa respectively. A compar- 
ison between the age of the oldest stars 
and the one that we have obtained by us- 
ing the best fit values of parameters, show 
that WDGP passes the age constraint very 
well. On the other hand we investigated 
the effect of the WDGP model on the ge- 
ometrical parameters of the universe such 
as the transverse comoving distance (co- 
moving angular diameter distance), the co- 
moving volume element and the angular 
size. And finally, we studied the effec- 
tive cosmological dynamics of the model 
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Abstract 

We investigate observational constraints on the normal branch of the warped DGP 
braneworld cosmology by using observational data from Type la Supernovae (SNIa), 
Baryon Acoustic Oscillations (BAO), Cosmic Microwave Background (CMB) and Baryon 
Gas Mass Fraction of cluster of galaxies. The best fit values of model free parameters are: 
Vtm = 0.240^^30 and = 0.000+°°^^ at la confidence interval by using Gold sample 
SNIa+CMB shift parameter+BAO+Gas mass fraction of baryons in cluster of galaxies. 
The results for essence sample SNIa combined with CMB shift parameter, BAO and 
Baryon Gas mass fraction correspont to: = 0.220+|];?7[J and n^.^ = 0.000+°-°25 at 68.3% 
confidence interval. We determine the age of the universe by using these best fit values. 
We also study the effective cosmological dynamics on the brane via an effective equation of 
state parameter and the deceleration parameter to conclude that an effective phantom-like 
behavior arises in this scenario. 



PACS: 98.80.-k, 98.80.Es, 95.36. +x 
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1. Introduction 

The accelerated expansion of the universe 
supported by recent observational data [1] 
is one of the most important discoveries in 
the last decade for the cosmology commu- 
nity. Within the framework of the gen- 
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eral relativity, the acceleration could be 
associated with the so-called dark energy, 
whose theoretical nature and origin are still 
unknown for theorists. Cosmological con- 
stant or vacuum energy with an equation 
of state parameter u = —1, is the most 
popular candidate for dark energy but un- 
fortunately, it suffers from some serious 
problems such as fine-tuning and coinci- 
dence problems. Therefore, a number of 

November 15, 2011 



models containing dynamical dark energy 
have been proposed as the mechanism for 
late-time cosmic speed up [2]. Some of 
them are quintessence, k-essence, phantom 
scalar field, chaplygin gas models and so 
on. Another alternative approach to explain 
the late-time cosmic speed up is modifica- 
tion of the geometric sector of the Einstein 
field equations leading to modified gravity 
[3]. In the spirit of modified gravity pro- 
posal, the Dvali-Gabadadze-Porrati (DGP) 
braneworld scenario explains the late-time 
accelerated expansion in its self-accelerating 
branch without need to introduce a dark 
energy component on the brane. In this 
scenario our universe is a 3-brane embed- 
ded in a higher dimensional flat space-time 
(bulk). The late-time acceleration is driven 
on large scales by leakage of gravity from 
the brane into the bulk ( see for instance 
[4] and references therein) . On small scales, 
gravity is bound to the brane and general 
relativity is recovered to a good approxima- 
tion. Indeed, the DGP model has two types 
of solutions (corresponding to two possi- 
ble embedding of brane in the bulk ): the 
self-accelerating, (DGP^'^^) branch and the 
normal, {DGP^-^) branch. The (DGP^^^) 
branch however, suffers from some insta- 
bilities such as ghosts [5] and cannot de- 
scribe the early stages of the universe evolu- 
tion properly. There are some extension of 
the DGP setup that provide relatively wider 
parameter spaces with richer phenomenol- 
ogy. One of these models is the warped 
DGP braneworld (WDGP) [6] which is a 
unified model of Randall- Sundrum II (RSII) 
braneworld scenario [7] and the DGP setup. 
In the RSII braneworld model, gravity mod- 
ifies in early ( high energy) epoches of the 
universe evolution. The warped DGP sce- 
nario also gives a self-accelerating phase in 
the brane cosmology. It is important to note 



that, the self-accelerating branch gives an 
effective equation of state that never can be 
less than —1 (always non-phantom behav- 
ior and therefore no crossing of the phan- 
tom divide line). The other branch of this 
scenario is the normal branch {DGP^~^) 
which doesn't self-accelerate but requires 
dark energy or modification of the induced 
gravity on the brane to explain the late 
time acceleration. It is possible to real- 
ize phantom-like effects (without phantom 
matter) in this normal branch via screening 
of the brane cosmological constant at late 
time [8]. 

This paper is devoted to explore obser- 
vational status of the normal branch of the 
warped DGP braneworld cosmology. Ob- 
servational constraints in DGP model with 
and without tension is investigated in Ref. 
[9]. We impose constraints on the model pa- 
rameters by using the several recent obser- 
vations such as distance measurements from 
type la supernovae (SNIa) from the Gold 
[10] and Essence [11] surveys, the baryon 
acoustic oscillations (BAO) measurement 
from the large-scale correlation function of 
the Sloan Digital Sky Survey (SDSS) [12], 
the position of the first peak of the cos- 
mic microwave background (CMB) from 
WMAP7 [13] and the baryon gas mass frac- 
tion of cluster of galaxies [14]. The struc- 
ture of the paper is as follows: in section 
2, we introduce the model and its cosmo- 
logical implications. In sec. 3 we explore 
the effect of the WDGP model on the ge- 
ometrical parameters of the universe. Sec- 
tion 4 includes the observational constraints 
on the model parameters space while in sec- 
tion 5 the detailed results are presented. In 
section 5 we perform a detailed comparison 
between age of the oldest objects of the uni- 
verse and the age result that obtained from 
the best fit values of our model parameter 



space. We also study the effective dynamics 
of the model in this section. Finally section 
6 is devoted to the concluding remarks. 



2. The Model 

The action of 
brane world model 
follows [6,15] 



the warped DGP 
can be written 



5 — Sbulk + 5; 



^brane j 



bulk 



2kI 



as 



(1) 



+ 



brane 



± 



brane 



(2) 

Here Sbuik is the action of the bulk, Sbrane 
is the action of the brane and S is the total 
action. with A = 0,1, 2, 3, 5 are coordi- 
nates in the bulk, while x'^ with /i = 0, 1, 2, 3 
are induced coordinates on the brane. k| is 
5-dimensional gravitational constant. ^^^R 
and ^^^Cjri are 5-dimensional Ricci scalar 
and matter Lagrangian respectively. 
is trace of the extrinsic curvature on either 
sides of the brane. Cbraneidapyip) is the ef- 
fective 4-dimensional Lagrangian. The ac- 
tion S is actually a combination of the 
Randall- Sundrum II and the DGP model. 
In other words, an induced curvature term 
is appeared on the brane in the Randall- 
Sundrum II model. Now we consider the 
brane Lagrangian as follows 

C-braneidap-, = — -R — A + L^, (3) 

where /i is a mass parameter, R is the Ricci 
scalar of the brane, A is the tension of the 
brane and is the Lagrangian of the other 
matter fields localized on the brane. We as- 
sume that bulk contains only a negative cos- 
mological constant, A5. With these choices. 



action (1) gives either a generalized DGP 
or a generalized RS II model: it gives DGP 
model if A = and A5 = 0, and gives RS II 
model if = 0. The generalized Friedmann 
equation on the brane is as follows [6] 



3/i2 



p + pQ{l + eA{p,a] 



(4) 



where e = ±1 is corresponding to two possi- 
ble branches of solutions (two different em- 
bedding of the brane) in this warped DGP 

.1/2 



model and A 



where A< 







1 - 2r] 
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V 



with < 
6^. By definition 



< 1 and 



mx 



Po 

A^/^ and 



6mg 
4 



m\ + 
ms = 



k^'^^^ . Also, Sq is an integration constant 
and corresponding term in the generalized 
Friedmann equation is called dark radiation 
term. We neglect dark radiation term in 
what follows. In this case, the generalized 
Friedmann equation (4) takes the following 
form 



a'' 



1 

V L 



,2 2r/p\i/2' 

p + po + £po[Ao^ 

Po ^ 



(5) 

where p = pm is the energy density of dark 
matter on the brane. 

2.1. Cosmological Implications 

In this section we study cosmological dy- 
namics on the DGP brane embedded in a 
warped bulk Manifold. To this end, we as- 
sume a flat FRW universe on the warped 
DGP brane. In this setup we can rewrite 
equation (5) as follows 




where 



\ is the DGP crossover scale. 



2m^ 

In the distance scale lower than this scale, 
gravity behaves as usual general relativis- 
tic one but in the distance scales higher 
than the crossover scale, gravity leaks to 
the extra dimension and this leakage leads 
to weakness of gravity in the large scales, so 
the universe expansion accelerates. The up- 
per sign in equation (6) corresponds to the 
self-accelerating branch of the model. Tak- 
ing the lower sign of this equation results a 
very interesting feature. Indeed, if we as- 
sume a model universe with standard cold 
dark matter (SCDM) with pm = Po„ (^) 
the accelerating behavior of the model can 
be recovered by rewriting Friedmann equa- 
tion (6) as follows 



+ A, 



where Ag// mimics the role of an effective 
cosmological constant on the brane (note 
that it is not actually a constant!) and it 
can be decomposed into two parts as follows 



A, 



A 
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A 



A5 
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The first two terms appeared in parenthe- 
sis on the right hand side of this relation 
could be considered collectively as a cosmo- 
logical constant term on the brane and the 
last term on the right hand side (the square 
root) screens the effect of the brane cos- 
mological constant in the same way as has 
been pointed out by Lue and Starkman [8]. 
In this situation, the effective cosmological 
constant A^ff on the brane increases with 
time due to dynamical screening effect, that 
is, reduction of the second term on the right 



hand side of (8) with cosmic time. In fact 
the normal branch of the model has the key 
property that brane is extrinsically curved 
so that shortcuts through the bulk allow 
gravity to screen the effects of the brane 
energy-momentum contents at Hubble pa- 



rameters H 



The screening effect is 



a result of leakage of gravity to the extra 
dimension at late times. 
For future purposes, it is useful to express 
the Friedmann equation (6) in a dimension- 
less form as follows 



H\z) 
HI 



where VLm = 



[1 + 2)3 + nx + fir, + fi 



As 



(9) 



POm. 
2 IJ2 
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3H, 



4f and Hn 



lOOh km/s/Mpc. Taking z = imposes 
a constraints on the model parameters as 
follows 



"(10) 

The general relativistic limit can be recov- 
ered if we set = (or = 0). In this 
case equations (10) implies that Q^ + ^x = 
1. Considering the tension of the brane 
cosmological constant this case leads to the 
ACDM cosmology. 

3. The effect of WDGP on the Ge- 
ometrical Parameters of the Uni- 
verse 

The cosmological observations are mainly 
dependent on the background geometry ( 
especially background spatial curvature) of 
the universe. So, in this section we study 
the effect of the WDGP model on the 
geometrical parameters of the universe. 



1. Comoving radial Distance 

One of the basic parameters in cosmology 
is the comoving radial distance. For an ob- 
ject with redshift 2; in a FRW background, 
this parameter can be expressed as follows 



Ha VP 



dz' 



H{z')/Ho 



(11) 

where J-" = {x, sinx, sinhx) for K = 
(0, 1, —1) respectively, i^' marks curvature 
of the spatial geometry and = '^jjjp- 
Figure 1 shows the radial comoving distance 
versus the redshift for different values of Vl^^ 
in a flat background. Clearly, increasing the 
values of Vtr^ results in a longer comoving 
distance. The mentioned quantity is a use- 
ful quantity in the analysis of the luminosity 
distances of Supernova type la. 
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Figure 1: Comoving radial distance versus the 
redshift for different values of & FRW 

universe. The unit of vertical axis is c/Hq. 

2. Angular size 

The apparent angular size of an object 
located at the cosmological distance is an- 
other important parameter that can be af- 



fected by the cosmological model during the 
history of the universe. If we take the ob- 
ject to lie perpendicular to the line of sight 
and to have physical extent the apparent 
angular size 9 is given by 



6 



D 



dAiz) 



:i2i 



where dA^z) = r{z)/{l + z) is the angular 
diameter distance which is a measure of how 
large objects appear to be in the universe. 
A key application of equation (12) is in the 
study of features of the cosmic microwave 
background radiation. The variation of ap- 
parent angular size AO in terms of the Az 
is given by 



% = H{z)r{z) 



(13) 



This relation is the so-called Alcock- 
Paczynski test. The advantage of the 
Alcock-Paczynski test is that in this case, 
instead of using a standard candle, we 
use a standard ruler such as the baryonic 
acoustic oscillation. Figure 2 shows Az/A6 
for different values of in a flat FRW 
background. As the figure shows, increas- 
ing of results in larger values of Az/ AO. 

3. Comoving Volume Element 

The comoving volume element is another 
geometrical parameter which is used in 
number-count tests such as lensed quasars, 
galaxies, or clusters of galaxies. The co- 
moving volume element can be expressed in 
terms of comoving distance and Hubble pa- 
rameter as follows 



/ 



dV _ r'^{z) 
dzdVt ~ H(z) 



(14) 



Figure 3 shows the comoving volume ele- 
ment versus the redshift for different values 



Figure 2: lS.z/ /^O versus the redshift for differ- 
ent values of in a flat FRW background. 

of VLr^ in a flat FRW background. This fig- 
ure indicates that the value of the comoving 
volume element increases with the increas- 
ing of fire- We note the quantities displayed 
in figures 1-3 are specified just by giving 
the value of fi^^ alone. This is because we 
are interested in the DGP character of the 
model. Nevertheless, in plotting these fig- 
ures we have used the values VLm = 0.24, 
Qf^^ = 1 and the other quantity Q\ is ob- 
tained via constraint equation, (10). 

4. Observational Constraints 

In this section, we constrain the model 
parameters of the warped DGP scenario by 
using the most recent observational data 
including the SNIa data measurements as 
given by the gold and essence samples, 
combined with the information from the 
BAO measurements by SDSS, the CMB 
shift parameter from WMAP7 observations 
and the baryon gas mass fraction. 

A: SNIa data 



Figure 3: Comoving volume element versus the 
redshift for different values of in ^ Acit FRW 
background. 

Here we use some new sample supernova 
type la data such as the gold sample com- 
piled in [10] and essence sample compiled in 
[11] to constrain the free parameters of the 
model. This observation directly measures 
the apparent magnitude m of a supernova 
versus its redshift z. We note that deference 
between SNIa samples is in their systematic 
errors and the range of redshift which the 
apparent magnitude is determined . 

The theoretical distance modulus is de- 
fined as 

fi^m-M ^5 logio dL + 5 logio (j~^^f^ + 

(15) 

where M is the absolute magnitude that is 
believed to be constant for all Type la su- 
pernovae and di is the luminosity distance 
that in general case can be expressed as fol- 
lows 

y/Wli] \ Jo E{z'-{ei})J 

(16) 

where {6i\ denotes the model parameters 
and E{z) is given by equation (9). We es- 
timate the best fit of the set of parameters 



{9i} by using a statistics, with 



N 

E 



Pobs{zj;{Oi}) - nthizj] {(^i})] 



(17) 

In this relation, is the number of SNIa 
data points which is different in several data 
sets, ^obs is the observed distance modulus 
and the aj is the uncertainty in the observed 
distance modulus, which is assumed to be 
Gaussian and uncorrelated so that the like- 
lihood is proportional to exp(— x^/2). The 
parameter M = Slogj^Q d ^^^jpc ) + 25 is a 
nuisance parameter and it is independent of 
the data and the data sets. Following the 
techniques described in Ref. [16], the mini- 
mization with respect to M can be made by 
expanding the Xsn equation (17) with re- 
spect to M as 

xlNiiOi}) = A- 2MB + M^C (18) 
where 

[Mobs(2j; {e^}) - {ej, A/ = o)]2 



Am}) = E ■ 



and 



Mobs (zj ; {f i }) - Mth (zj ; {0i }, M = 0) 



j=l 



C ■- 



^ 1 

3 = 1 J 



(20) 



Equation (18) has a minimum for M = 

§ at xiAm = Am) - 

Using this equation the best fit values of 
model parameters as the values that mini- 
mize XsNiaii^i}) can be obtained. 

For the likelihood analysis we marginal- 
ize the likelihood function L ~ exp(— x^/2) 
over h. We adopted Gaussian priors such 
that h = 0.705 from the WMAP7 [13]. 
Table 1 summarizes these priors. 



Table 1: Priors on the parameter space used in the 
Ukehhood analysis. 



Parameter 


Prior 




0.000 - 1.000 Top Hat 




0.000 - 1.000 Top Hat 




-1.000 - 1.000 Top Hat 


h 





(19) 



B: CMB shift parameter 

We use the CMB data from WMAP7 ob- 
servation that includes the shift parameter 
TZ and the redshift of the decoupling epoch 
. The shift parameter TZ relates the angu- 
lar diameter distance to the last scattering 
surface, the comoving size of the sound hori- 
zon at = 1091.3 and the angular scale of 
the first acoustic peak in the CMB power 
spectrum of the temperature fluctuations. 
The CMB shift parameter is approximated 
by (see for more details Ref. [16]) 



7^ 




(21) 



where dA^z) is the angular diameter dis- 
tance defined by equation (12). The con- 
straints on a typical model using CMB shift 
is obtained from minimization of the quan- 
tity 



XCMB 



(22) 

'CMB 

where TZots = 1-725 is the observed value of 
the CMB shift parameter performed from 
WMAP7 observation and its corresponding 
la error is <Jcmb = 0.018 [13]. Also TZth 
corresponds to the theoretical value of shift 
parameter calculated from equation (21). 

C: BAO observation 
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The baryonic acoustic oscillation (BAO) 
peak detected in the SDSS luminous red 
Galaxies (LRG) is another tool to test the 
model against observational data. BAO are 
described in terms of a dimensionless pa- 
rameter 



■Aizsds3;{di}) = vf^ 



Hodl{zsdss;{9i}) 



H{Zsdsa; {^i})2sdss(l + Zsdas) 



1/3 



(23) 



The fo^^ the BAO is given by 



a 



(24) 



sdss 



•^obs 0.469 ( Q gg 



^obs 
-0.35 



The observed value from the LRG is 

± 0.017 measured 

at Zsdss = 0.35 [20]. Here Ug = 0.963 is 
the spectral index as measured by WMAP 
seven years observations [13]. 

D: Gas mass fraction of cluster of 
galaxies 

Another cosmological test to constrain 
the parameters of the model arises from 
baryon gas mass fraction of cluster of galax- 
ies for a range of redshifts 



M, 



gas 



gas 



Mi 



(25) 



tot 



The basic assumption underlying this 
method is that the baryon gas mass frac- 
tion in clusters is constant, independent of 
redshift. This method can give a constraint 
to the geometry of the universe with the re- 

3 

lation Sgas oc d\ under the assumption that 
this fraction should be approximately con- 
stant with redshift. Following [17] (see also 
[16]), the expression for gas mass fraction 
is given by 



lsfaU^j;{e,})~sl'u^j;{e,})]^ 



(26) 



where Sgas is a dimensionless parameter de- 
fined as 



S, 



b Qh ( dJl' 



gas 



1 + l3D.rn\ dA{z 



lat I 



(27) 



In this relation dJ)^^ is the angular diame- 
ter distance to a cluster in the test model 
which is assumed to be SCDM (cold dark 
matter) in this case, and 6 is a bias fac- 
tor suggesting that the baryon fraction in 
clusters is slightly lower than for the uni- 
verse as a whole. Also 1 + /3 is a factor 
taking into account the fact that the total 
baryonic mass in clusters consists of both 
X-ray gas and optically luminous baryonic 
mass (stars), the latter being proportional 
to the former with proportionality constant 
/3 ^ 0.19^/5^ [16-18]. The nuisance param- 
eter f = TTaTT^ should bc marginalized via 
expanding the x^as equation (26) with re- 
spect to ^ which gives 



X 



gas 



K - 



Y 



(28) 



where 



K 



<7^ 



(29) 



and 



W 



cobs I 



z,;{9,}).Sll{z,;{9,},i = l) 



Y 



^2 



(30) 



We use the 26 cluster data [14] for 5°^^ to 
obtain the best fit parameters of the model. 

It is important to note that the above 
observational data are uncorrelated since 
they are given by different experiments and 
methods. So, we can construct a joint anal- 
ysis as 



Xt 



ot 



XSN + XCMB + XSDSS + Xlas ■ (31) 



5. Results 



With these prehminaries, we have ob- 
tained the best fit parameters of the normal 
branch of the WDGP model for SNla data 
(gold and essence datasets), the joint anal- 
ysis of the SNIa and CMB, the combined 
analysis of the SNIa, CMB and SDSS and 
finally the joint analysis of the total datasets 
including cluster galaxies gas mass fraction. 
Table 2 shows the results of the observa- 
tional constraints on the free parameters of 
this model. Figure 4 shows the marginal- 
ized relative likelihood with respect to pa- 
rameter and fitted with SNIa gold 
sample, SNIa+CM]3, SNIa+CMB+SDSS 
experiments and SNIa+CMB+SDSS-|-gas 
mass fraction observations. In figure 5, we 
repeat these stages for Essence sample of 
SNIa experiments together with other data 
sets. We plot contour maps of ver- 
sus Qm in figure 6 for Gold Sample and 
combined observational data sets. In fig- 
ure 7 we plot corresponding contour maps 
for Essence Sample and combined observa- 
tional data sets. 
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Table 2: The fitting results for WDGP model by using the SNIa (Gold and Essence Samples), SNIa+CMB 
and SNIa+CMB+SDSS and SNIa+CMB+SDSS+glaxies clusters gas mass fraction experiments in a flat 
background. 



Observation 










Xiiin/^d.o.f 


Agc(Gyr) 


SNIa( Gold Sample) 




u.iou_o.i5o 




0.009 


0.923 


64.97±^^:^^1! 


SNIa(Gold)+CMB 


0.280±U:\'^lj 




0.631±U:f/^, 


0.939 


0.943 


14.38+-- 


SNIa(Gold)+CMB+SDSS 


260+" '^'^" 

U.ZOU_o.l20 


0.000+°°i2 




0.839 


0.992 


i4.5i±^:ti 


SNIa(Gold)+CMB+SDSS+GAS 


U.Z'±U_o.l30 


000+0014 




0.999 


0.992 


14 «f;+i+t)i 

i4.00_3 942 


SNIa( Essence Sample) 


U.UZU_o.015 


170+0.020 


lfi1+0-0V4 


0.009 


1.032 


78 QO+56.018 
(0.00_42 390 


SNIa(Essence)+CMB 


2-50+"-"^" 

U.ZOU_Q J^gQ 


000+0019 


n fiQI +O il'i' 


0.539 


1.030 


14.86+--^ 


SNIa(Essence)+CMB+SDSS 


220+" "^'^ 


00 +0 022 


fiq7+".i4b 
u.Dy(_o.095 


0.689 


1.054 


15.071-- 


SNIa(Essence)+CMB+SDSS+GAS 


220+'^ "^" 

^■^^"-0.170 


0.000 +0 025 


710+0.187 


0.119 


1.044 


14.96±^:i^ 




Figure 4: Marginalized relative likelihood with respect to parameter (left) and fi^c 

(right) fitted with SNIa Gold Sample, SNIa (Gold)+CMB, SNIa(Gold)+CMB+SDSS and 
SNIa(Gold)+CMB+SDSS+ galaxies clusters gas mass fraction experiments. The intersection of 
each curve with horizontal solid line is corresponding to the bound with la confidence level. 
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Figure 5: Marginalized relative likelihood with respect to parameter ilm (left) and ^Ir^ (right) 
fitted with SNIa Essence sample, SNIa (Essence)+CMB, SNIa (Essence)+CMB+SDSS and SNIa 
(Essence)+CMB+SDSS+galaxies clusters gas mass fraction experiments. The intersection of each curve 
with horizontal solid line is corresponding to the bound with Icr confidence level. 




0.1 0.2 0.3 0.4 0.1 0.2 0.3 0.4 0.1 0.2 0.3 0.4 0.1 0.2 0.3 0.4 



Q Q Q. £1 

m m m m 

Figure 6: Contour maps of ^Ir^ versus for WDGP model with la (68.3%) and 2a (95.5%) confi- 
dence levels using SNIa (a), SNIa+CMB (b), SNIa+CMB+SDSS (c) and SNIa+CMB+SDSS+Gas (d). We 
used the Gold Sample for SNIa experiment. 



11 



0.2 



0.15 



a" 0-1 



0.05 






(a) 






V\ 




\ 





0.05 



0.04 



0.03 



0.02 



0.01 




0.05 



0.04 



0.03 



0.02 



0.01 




0.06 



0.04 



0.02 




0.1 0.2 0.3 0.4 0.1 0.2 0.3 0.4 0.1 0.2 0.3 0.4 0.1 0.2 0.3 0.4 
Q. CI CI Q. 



Figure 7: Contour maps of ilrc versus 0,m for normal branch of WDGP model with la (68.3%) 
and 2a (95.5%) confidence levels using SNIa (a), SNIa+CMB (b), SNIa+CMB+SDSS (c) and 
SNIa+CMB+SDSS+Gas (d). We used the Essence Sample for SNIa experiment.. 



12 



5.1. The Age of the universe 

The age of the universe for an expanding 
universe in a flat background is given by 

to= I dt= I (32) 



.0 + 

The WMAP collaboration [13], using the 
AC DM model, quotes to = 13.75 ±0.11 
Gyr. In our model H{z) is given by Eq. 
(9). A model- independent limit of to can be 
obtained from the age of the oldest stars, 
tstar- Since the age of the universe should be 
greater than the age of the oldest stars, this 
limit can be used as another probe to con- 
strain our model ( see for instance [19,20]). 
Studies on the old stars [21] suggest an age 
of 13^^ Gyr for the universe. Richer et. 
al. [22] and Hansen et. al. [23] also pro- 
posed an age of 12.7 ± 0.7 Gyr, using the 
white dwarf cooling sequence method. Re- 
cently, Frebel et al. [24] have reported the 
discovery of HE 1523 — 0901, with an age 
tstar = 13.4±0.8(l(T)±1.8(s?/st) Gyr. With 
these statistical and systematic errors, we 
adopt tstar > 12 Gyr for constraining our 
model [21]. Table 2, shows that the age 
of the universe is always greater than tstar 
for all used data sets. However, the SNIa 
data sets (Gold or Essence Samples) alone 
lead to very great to which is not reliable 
and it seems that the pure SNIa test cannot 
give a good constraint on the model param- 
eters space. This feature is evident also in 
our previous relative likelihood figures and 
has its origin in the geometric nature of the 
SNIa data sets. On the other hand the com- 
bined analysis of SNIa with other data sets 
improves the constraints and gives a more 
reasonable results. 

5.2. Effective Dynamics 

It has been shown in our previous work 
[15] that the Warped DGP scenario has a 



phantom-like behavior without need to in- 
troduce a phantom matter neither in the 
bulk nor on the brane. To investigate the 
effective dynamics of the model we study 
the effective density defined by 



Peff = A4 



2rl 



1 + 4r2(^ 



3/i2 



P + A As 
3 
(33) 

and an effective equation of state parameter 



-1 



Peff 



(34) 



So, the effective equation of state of dark 
energy can be expressed as follows 



= -1- 



(35) 

Using the best fit parameters obtained 
here and summarized in table 2, in figure 
8 we plot ujeff versus the redshift for com- 
bined SNIa+CMB+SDSS+galaxies clusters 
gas mass fraction data set for Gold and 
Essence Samples. 

As the figure shows, in low redshifts and 
especially at present, dark energy exhibits a 
phantom-like acceleration with Ueff < — 1 ( 
it is easy to show also that in this model 
H < and Pe// > 0, see [15] for details). 
In future { z < 0) the effective equation of 
state parameter approaches the cosmologi- 
cal constant line Ugff = — 1. It is important 
to note that although Wg// has a phantom- 
like dynamics in this warped DGP scenario, 
ojtot doesn't show such a behavior. The ex- 
pression for Wtot is given by 



^tot 



(36) 
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Figure 8: The effective equation of state 
parameter versus the redshift fitted with 
SNIa+CMB+SDSS+galaxies clusters gas mass 
fraction with Gold and Essence Samples re- 
spectively. 
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Figure 9: Total equation of state pa- 
rameter versus the redshift fitted with 
SNIa+CMB+SDSS-h galaxies clusters Gas 
mass fraction with Gold and Essence Samples 
respectively. 



Indeed, the total equation of state param- 
eter remains Quintessence-like [wtot > —1) 
as figure 9 shows. 

Another dynamical quantity of interest 
is the deceleration parameter which is ex- 
pressed as g = 



HoE^ 



I where 



E 



1 - 



(37) 



This equation implies that for all values 
of 2;, E < 0. As an important result, the 
deceleration parameter could be never less 
than —1. Consequently, there is no super- 
acceleration or big rip singularity in this 
model. Figure 10 shows the variation of the 
deceleration parameter versus the redshift 
with two combined data sets. It is clear that 
the deceleration parameter reduces by red- 
shift towards recent epoch and in the future 



it approaches q = —1. 

6. Summary and Conclusion 

The normal branch of the warped DGP 
brane scenario has several fascinating prop- 
erties. Its cosmological dynamics does not 
suffer from instabilities such as ghost, it can 
exhibit a phantom mimicry and it leads to 
the famous ACDM scenario in certain lim- 
its. In this regard, we have studied cosmo- 
logical constraints imposed on this model 
from observational data such as type la su- 
pernova data from the Gold and Essence 
surveys, the baryon acoustic oscillations 
(BAO) measurement from the Sloan Dig- 
ital Sky Survey (SDSS), the cosmic mi- 
crowave background (CMB) and the baryon 
gas mass fraction of clusters of galaxies. 
We find that the best fit values of model 
free parameters are constrained to: fi^ = 
0.2401°;°^° and fi,^ = 0.000+° °i^ at la con- 
fidence interval by using Gold sample SNIa 
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0.000+°-°2^ at 68% confidence in- 



Figure 10: Variation of the deceleration param- 
eter versus the redshift. The parameters values 
used to produce this figure are taken from table 
2. 



+ CMB shift parameter + BAO + gas mass 
fraction. The same analysis just including 
essence sample SNIa instead of Gold sam- 
ple correspond to: 



0.2201°:??° and 



terval. Furthermore, using these best fit 
parameters, we obtained the age of the 
universe U.SbtHH and UMtHH at la 
confidence interval by using SNIa+CMB 
shift parameter+BAO+gas mass fraction 
for gold sample and essence sample of SNIa 
respectively. A comparison between the age 
of the oldest stars and the one that we have 
obtained by using the best fit values of pa- 
rameters, show that WDGP passes the age 
constraint very well. On the other hand we 
investigated the effect of the WDGP model 
on the geometrical parameters of the uni- 
verse such as the radial comoving distance, 
the comoving volume element and the an- 
gular size. And finally, we studied the ef- 
fective cosmological dynamics of the model 
via effective equation of state parameter 



and the deceleration parameter by using the 
best fit parameters of the combined analy- 
sis of SNIa+CMB+SDSS+galaxies clusters 
gas mass fraction with Gold and Essence 
Samples. This analysis confirms that the 
WDGP model has a phantom-like behavior 
without need to introduce phantom fields 
neither in the bulk nor on the brane. It is 
important to note that in this scenario the 
effect of the bulk's warped geometry (en- 
coding in A5) is that by increasing the abso- 
lute value of the bulk cosmological constant 
the phantom-like behavior reduces. In other 
words, incorporation of A5 leads to a reduc- 
tion of the effective phantom nature of the 
model in comparison with pure DGP case ( 
see Ref. [15]). 
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